DESY 05-154 



Hot QCD, /c-strings and the adjoint monopole gas model 



Chris P. Korthals Altes 
Centre Physique Theorique an CNRS 
Case 907, Campus de Luminy, F13288, Marseille, France 

Harvey B. Meyer 1 
DESY 
Platanenallee 6 
D-15738 Zeuthen 



Abstract 

When the magnetic sector of hot QCD, 3D SU(iV) Yang-Mills theory, is described as a dilute 
gas of non-Abelian monopoles in the adjoint representation of the magnetic group, Wilson 
loops of A/"-ality k are known to obey a periodic k(N — k) law. Lattice simulations have 
confirmed this prediction to a few percent for N = A and 6. We describe in detail how the 
magnetic flux of the monopoles produces different area laws for spatial Wilson fc-loops. A 
simple physical argument is presented, why the predicted and observed Casimir scaling is 
allowed in the large-iV limit by usual power-counting arguments. The same scaling is also 
known to hold in two-loop perturbation theory for the spatial 't Hooft loop, which measures 
the electric flux. We then present new lattice data for 3D N = 8 fc-strings as long as 3'fm' 
that provide further confirmation. Finally we suggest new tests in theories with spontaneous 
breaking and in SO (An + 2) gauge groups. 
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1 Introduction 



The title of this paper may sound to most practitioners of lattice gauge theory and hot QCD of 
a somewhat esoteric nature. And on the other hand aficionados of the beauty of non-Abelian 
monopoles |S1 IH Ej [HI E] pH CHI ma y reflect on the title as being heretic, since non-Abelian 
monopoles have so far withstood the traditional approach that has been implemented success- 
fully for 't Hooft-Polyakov monopoles: as yet, nobody has come up with a viable construction 
of a classical solution, that is then quantized by semi-classical methods. Only recently a con- 
struction of non-Abelian fluxes in a low energy field theory version has been accomplished H . 
These are models that relegate the intricacies of non-Abelian monopoles to their high energy 
sector, and manage to construct explicit non-Abelian fluxes in the low energy sector. 

Somewhat analogously, we will forget about the intricate nature of individual non-Abelian 
monopoles and assume that a gas of such objects has relatively straightforward properties. 
That will allow us to compute and interpret in a simple-minded way the average behaviour of 
magnetic flux loops, that is, spatial Wilson loops [TH] . Such spatial loops have been measured 
in lattice simulations by Teper's group [231211 in a wide temperature range, thus showing that 
the predictions of the model can be tested from first principles. 

At temperatures well above the critical T c , the temporal extent of the system becomes 
negligible, and we are left with a three-dimensional system. This implies that the tension of 
the spatial loop at such very high T also bears the interpretation of a three-dimensional string 
tension, due to a chromo-electric flux tube. In other words, our model is also indirectly a model 
for confinement in 2+1 dimensional gauge theories. 

Generally speaking, in non-Abelian SU(iV) gauge theories in three and four dimensions, 
the chromo-electric flux between two static colour sources arranges itself so as to produce a 
linearly rising potential. This naturally suggests a flux-tube configuration and leads to the 
string picture of confinement. While in SU(3) there is only one string tension, that of the 
string appearing between charges in the fundamental representation, in SU(iV > 4) there are 
[N/2] independent stable '/c-strings' which are protected from screening by the center-symmetry 
Z(N). 

The picture that we propose for the origin of the area laws of the spatial Wilson /c-loops, and 
hence for 3d fc-strings, is rooted (perhaps paradoxically) in high temperature 3+1 dimensional 
QCD and involves a gas of screened non-Abelian monopoles - or rather "magnetic quasi- 
particles". We prefer the latter terminology, since it stresses that our monopoles need not be 
eigenstates of the Hamiltonian but are rather collective modes of the plasma. The objects that 
we shall describe in the 3d gauge theory are the dimensionally reduced versions of these modes, 
much in the same way as Polyakov's 'pseudoparticles' flj in the 3d Georgi-Glashow model 
are the descendants of the t 'Hooft-Polyakov monopoles |2] living in the 4d Georgi-Glashow 
model. The non-Abelian Stokes theorem |H] establishes a connection between spatial Wilson 
loops and the magnetic flux in the plasma; which in our model is induced by the magnetic 
quasi-particles. That is, schematically, how we are able to make predictions for 3d /c-string 
tensions. 

Of course, fc-strings are also interesting in their own right. Since they are perfectly stable, 
their tension ratios can be used to discriminate unambiguously between models of confinement. 
In what follows, without giving a comprehensive view of the latter, we put our model in 
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perspective with respect to a broader class of such models. 

Our adjoint monopole gas model jTHl EE] is related to the dual-superconductor picture of 
confinement [I]. The latter would naturally predict the presence of monopoles in the plasma, 
as manifestations of the condensate at low T. It is a natural generalization of the seminal idea 
of 't Hooft ^B], that Abelian monopoles Bose-condense in the ground state, and are transient 
states in that they won't show up in the spectrum of the Hamiltonian. In the hot deconfined 
phase they should populate the ground state, just like gluons. To explain the /c-loop tensions 
in the hot phase is however non-trivial because the number of different species of Abelian 
monopoles is too small (N — 1 for SU(iV)). 

There is the elegant caloron solution to the equations of motion [23J. It is a periodic 
instanton with a Higgs-like background furnished by the non-trivial value of the Polyakov 
loop. This gives rise to N monopoles (a fundamental multiplet) inside the caloron. Could 
these be related to the quasi-particles that we are invoking? It may be [23] that a ^ high 
enough temperatures the monopoles inside an individual caloron start to "deconfine" and are 
able to move freely from one to another caloron, much in the same vein that gluons can freely 
move from one glueball to another at high T. However free monopoles in the fundamental 
multiplet can not explain the observed Casimir scaling [21]. Nevertheless, as explained in the 
next section, even at asymptotic temperatures we are actually facing strong coupling when we 
try to explain the spatial Wilson loop behaviour. It could well be that this strong coupling 
favours binding into adjoint monopoles (while binding into singlets is statistically disfavoured 
at large N). Non- Abelian monopoles in the adjoint representation furnish precisely the correct 
number of species to explain the observed Casimir scaling, as shown in earlier work 19 and 
in section [Ol below. 

The ratios of A;-string tensions are also tests for formulations of SU(iV) gauge theories 
derived from fundamental string theory. Examples of the latter are the MQCD framework jUl] 
and the AdS/CFT calculations in Ref. [12] for D = 3. The latter give Casimir scaling for N 
large and k of order N; our model predicts Casimir scaling for any value of N. 

The MQCD framework gives a sin(/c7r/iV) law for the k-tension, implying in particular that 
the tension ratios o~k/o~\ have 1/iV 2 corrections. An elegant paper by Gliozzi [22] provides a 
simple geometric interpretation for the sine law. He shows that in the cold phase the sine law 
is the borderline between formation of Z(N) symmetric static baryons (no k > 2 flux tubes 
involved) and formation of static baryons with k = 2 or higher flux tubes (it is assumed that 
arbitrary short flux tubes have the same tension as long flux tubes). 

Casimir scaling and the sine law both predict that o~\.fo~\ — ► at large N, fixed k; in other 
words, a fc-string is a collection of k non-interacting fundamental strings in the planar limit 
N — > oo. Casimir scaling however attributes a binding energy to these k strings of order 
1/N, while if the sine law is correct, this energy is only 0(1/N 2 ). Recently there has been 
a discussion 39 on the conflict of 1/N corrections with standard 1/N power-counting rules, 
based on the assumption that all representations with a given A/"-ality k produce the same 
tension. We point out in section [B] that this analysis neglects mixing effects between reducible 
representations which are of order 1/N and which lower the energy of the lightest string by an 
amount of that order. Earlier work on strong coupling expansions [13] corroborates our general 
argument. More recently, analytic calculations of the tension for 't Hooft loops [201 EI] have 
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been shown to lead to the same k(N — k) scaling law. 

As already mentioned, lattice calculations have been carried out [221 I2H1 I2H] i n three and 
four dimensional SU(iV) gauge theory to determine the ratios of the A;-string tensions to the 
fundamental string tension. Here we study the /c-strings in 3d SU(iV) gauge theories, present- 
ing new data for their tension ratios obtained for the gauge group SU(8) and combining the 
new information with previously obtained SU(4) and SU(6) data |2S]- The numerical advan- 
tage of searching for the effect of monopoles on Wilson loops at high T is that the relevant 
simulations are three-dimensional; needless to say, to obtain the same accuracy, the amount of 
computational effort is considerably lower for the 3d simulations employing the reduced action. 

By the same token, given such an accuracy for the 3d lattice data, it is useful to know to 
what accuracy in the coupling g(T) the dimensionally reduced actions reproduce the full 4d 
QCD result. For the case of three colours one knows [SB] that the 3d results for the string 
tension reproduce the 4d lattice data up to 1.1T C through the running of g(T) up to and 
including two loops. 

The lay-out of the paper is as follows. We start with section El on how the problem of the 
residual strong interactions in hot QCD is attacked quantitatively - by dimensional reduction. 
In section El we review briefly non-Abelian monopoles. We then derive in section [Ol a Stokes 
type formula for the spatial Wilson loops that permits us to quantify the effect of the putative 
non-Abelian monopoles in section 14.41 Then follows section El on strings in higher represen- 
tations where our arguments on the 1/N counting are exposed, and the lattice calculation is 
presented in section El Finally we compile and discuss the lattice data accumulated so far 
(section E31) and the paper ends with a general conclusion (section EJ). 



2 High temperature QCD 

This section is meant to introduce the reader into the essentialia of hot QCD, and to motivate 
the model. 

At temperatures well above T c asymptotic freedom drives the running coupling g(T) down 
to zero. On the other hand the average density of gluons is the Bose-Einstein density ube{p/T) 
(jp = \p\ is the momentum of a gluon). Because of this density the coupling in the plasma has 
to describe stimulated emission and equals 

9 1 = 9 2 n B E{p/T) = g 2 — — -. (1) 

expp/T — 1 

This leads to a picture of a gluon plasma, where one has to distinguish three scales: 

• hard gluons with momentum p of order T, interacting weakly, g 2 st = 0(g 2 ). 

• soft gluons with momentum p of order gT, still interacting weakly g 2 t = O(g). 

• ultra-soft gluons with momentum p = 0(g 2 T), interacting strongly, g 2 t = 0(1). 

Thus, in spite of asymptotic freedom, there is a strongly interacting sector left. Strongly 
interacting because the large population of ultra-soft energy levels pushes the coupling up [12] ■ 
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Thus, at these length scales, semi-classical methods are unlikely to apply, as we argued in the 
introduction. 

The hard gluons are familiar from the Stefan-Boltzmann form for the pressure. The hard 
gluons cause Debye screening mo~ gT of the force between electric test charges. All this has 
been known for long from electrodynamic plasmas. 

A new feature becomes apparent for the non-Abelian plasma at scales g 2 T. It is the 
screening of the magnetic force (thm ~ g 2 T) between two static magnetic test charges. In 
electrodynamic plasmas no static magnetic screening exists. Magnetic screening not only 
occurs at arbitrary high temperatures, it persists at arbitrary low temperatures, where the 
electric screening has disappeared and has turned into electric confinement. It is a hallmark 
of the non-Abelian system, and hints at a magnetic activity for all temperatures ^2] • 

One purpose of this paper is to understand and test a specific model [Tj5] for the strongly 
interacting sector. We will state its assumptions at the end of this section. 

2.1 Dimensional reduction at high T 

In this section we give a fast review of how one computes equilibrium properties of the plasma 
in a systematic way. The problem of strong coupling at large distances is dealt with through a 
sequence of effective actions |3S] ■ It is the last and strongly interacting effective ( "magnetic" ) 
action that our monopole model approximates. 

By integrating out the hard modes in the QCD action one produces an effective 3d action 
called Seqcd- If one accepts to have an accuracy of 0(g(T) 4 ) this electrostatic action is the 
superrenormalizable action in terms of the static potentials. The form of our effective action 
Seqcd is dictated by all symmetries, global and local, of the original QCD action, which are 
respected by the integration process. That implies all the symmetries we knew already, except 
that the electric term in the static action will have no d$A term. So A Q appears as an adjoint 
Higgs term in our 3D gauge theory. The electrostatic QCD action density reads: 



Because of R- conjugation invariance (Aq — ► — Aq) the electrostatic action must be even in Ao. 
For SU(2) and SU(3) the second quartic term is identically zero. 

The parameters in this 3d action are the coupling g^, the electric mass tjie and the 4-point 
couplings A. All of them are expanded in powers of the QCD running coupling g 2 (T), and all 
of them are now known to 0(g A ) [^niEZI- The electric mass coincides with the Debye screening 
mass m D = 2-^-T 2 to one loop order. The 4-point couplings start with the fourth power of 
g(T). It is customary [23] to express all these parameters in terms of the dimensionful scale 
gE, x = \e/9ei similarly for x, and finally y = rn 2 E /g E . For large T the x variable becomes 
equal to g 2 , the xy variable approaches a constant. 

In the limit where the electric mass gT becomes very large compared to the coupling 

g 2 E = g 2 T one can integrate out this mass scale and obtain magnetic QCD: 



= Tr {(D(A)Ao) 2 } + m|Tr {A 2 } + A £ (Tr {A 2 }) 2 + 
+ X E [(Tr {A }) 4 - ~(Tr {A 2 }) 2 ] + ^Tr {**} + 8C E . 



(2) 



±Tt{F 2 } + 5£ m . 



(3) 
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The coupling parameter in this Lagrangian is called gu and it can be expressed in terms of 
the electric coupling g E [38] : 



9% 17 / g% 



) 



2 



2 2 i 

9m = 9e 1 



IQnrriE 512 \7rm£ 



(4) 



One should realize that in the pure N = 3 Yang-mills theory there is only one parameter: 
the coupling g{T / 'At)- This means there is a relation between x and y, where the physics of 
the plasma is: 



This action serves to compute the leading contribution to magnetic quantities like the 
spatial Wilson loop a, or the magnetic screening tum at very high T. For dimensional reasons 
both are proportional to g 2 M . The corrections are very small, and this seems to be a general 
feature of this type of corrections 08] • On the other hand the corrections due to hard modes 
in two loop approximation are appreciable [38]. It turns out that one can extrapolate the 3d 
result for the Wilson loop to about 1.1T C just by using this 2-loop running of the coupling 
gE, as a very good approximation to the 4d results. Quite likely the same is true for the 
magnetic screening length Im or magnetic screening mass ttim = Iaj, which is defined from 
the correlation of a heavy monopole pair: as for the spatial tension, its dominant contribution 
comes from the 3d magnetic sector. 

Summing up: computing magnetic quantities at T 3> T c , in 3d magnetostatic QCD, is 
sufficient to know them all over the deconfined phase by simply using the two loop running of 
the coupling. This means that some salient features of our model for the magnetic sector (see 
next subsection) are valid for all of the deconfined phase. 

2.2 Magnetic quasi-particle model for the magnetic sector 

The magnetic sector is governed by 3d Yang-Mills theory. For the physics over distances larger 
than the magnetic screening length we make three Ansaetze: 

1. The interaction for the magnetic gluons is so strong that they bind in lumps. 

2. The lumps are dilute. 

3. The lumps are non-Abelian monopoles. 

Their size is on the order of the magnetic screening mass tum = 0{g 2 T). And so is their 
inter-particle distance, or their density Um- The ratio of the two corresponding volumes is the 
diluteness 



As the coupling drops out in this ratio there is no parametric reason that the diluteness is 
small. 




for N = 3. 



(5) 



5 = n M /m 3 M . 



(6) 
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From these Ansaetze follows from simple dilute gas arguments (repeated in section 14.4)1 
that the tension a of the spatial Wilson loop equals: 



1 



(7) 



m M 



There is a group factor in front (which will be discussed in section l4"3)l . and the the corrections 
are in powers-not necessarily integer- of the diluteness. 

So the diluteness is known, once the tension and the magnetic screening mass are known 
from lattice measurements. Its smallness is a dynamical effect giving a value of about 0.05 
(with a correction of 0(1/N 2 ) j21j) in the hot phase for SU(iV) groups: it is given by the string 
tension in the 3D gauge theory in units of the lightest glueball mass. Clearly it is gratifying 
to have an - admittedly empirical - justification for the diluteness being small. 

It is instructive to compare our dilute gas of composite lumps with radius Im to the dilute 
gas one finds in the usual weak coupling plasmas. There the lumps are point-like particles, and 
the Debye screening length Ijj is large with respect to the inter-particle distance, i.e. l z D n 3> 1, 
the weak coupling plasma condition. For hard gluons one has n ~ T 3 and l^ 2 ~ g 2 and the 
plasma condition is fulfilled. 

We want to close this section with a brief comment. It is tempting to go down from infi- 
nite temperature to finite temperature, and consider the lumps as magnetic quasiparticles, or 
collective excitations of the plasma. From our theoretical knowledge of the magnetic screen- 
ing [21] we know that at T = the screening mass equals the lowest glueball mass in the 
4D gauge theory; and the spatial Wilson tension equals the string tension at T = 0. Lattice 
simulations j25| find the diluteness at zero temperature, as given by the string tension in units 
of the lightest glueball mass in the 4D gauge theory, is still small, on the order of 0.09 for 
N large! Thus our dilute gas stays dilute when lowering the temperature. At some temper- 
ature T q the Bose-Einstein statistics takes over (where the ratio of magnetic screening and 
de Broglie thermal wave length T _1 become on the same order). And, in the spirit of dual 
superconductivity, BE-condensation is then marking the transition to the confined phase. 

3 Non-Abelian monopoles 

The magnetic sector of hot QCD has magnetic lumps through the strong (g 2 ~ 1) binding of 
magnetic gluons. Very specifically we do not have a Higgs field at our disposal to define the 
U(l) field strengths. The question is whether other than 't Hooft-Polyakov [2] monopoles can 
be formed under such circumstances. The answer is not known to date. But if they are there 
they must obey a Dirac condition. 

In 1977, Englert, Windey, Goddard, Olive and Nuyts [Sj analysed precisely such hypothet- 
ical monopoles in an unbroken gauge theory, and formulated the generalized Dirac condition. 
This condition is the following. Let B be a matrix in the SU(iV) Lie-algebra. Let the colour 
magnetic field B be given far away from the monopole by: 



B 



(8) 



B = gr 



4 7rr 3 
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The Dirac condition then reads: 

expigB = l. (9) 

This condition has to be fulfilled for any matter field that couples to the gauge field. Obviously 
we can take B to be diagonal. Note also that for U(l) we get the expected result gB = 2nn. 

For any simple Lie group one has the orthogonal set of diagonal generators H = (Hi, . . . , H r ) , 
with r the dimension of the Cartan subalgebra. The remaining orthogonal generators are 
E a = EL a . The roots a = (ai, . . . ,a r ) are given by: 



1 






■E± a 


a 




a 






■ H. 


a 2 



[H, E a ] = aE a 

[E a , Ep] = (a + (3)E a+ p if a + (3 is a root 

[E a ,Ep] = otherwise 

[E a ,E_ a ] = a-H. (10) 

These definitions imply a common normalization Tr {E a E^ a } = Tr{Hf}. In physics we are 
used to have it equal to 1/2. 

We define now the coroots a — In terms of those the group admits a set of SU(2) 
subgroups (like the familiar I, U and V spin in SU(3)) for any root a, denoted by SU(2) a . One 
gets them by projecting H on the coroots and using the E a . More precisely: 

E± a 

H a = ^-H. (11) 

Crucial is now that the matrices H a , being homogeneous in the roots and the Hi, are inde- 
pendent of the normalization of the matrices Hi. Hence they have eigenvalues, which are pure 
numbers. What are those? 

The commutation relations that normalize H a follow from Eq. with the result: 

[E a , E_ a ] = H a . (12) 

The weights of an irreducible representation are given by the eigenvalues of the diagonal op- 
erators. So if the carrier vectors V(k) diagonalize the representation we can define the weights 
w k by : 

Hv( k ) = w k v {k) . (13) 

Now a theorem on Lie algebras JT] tells us that 2H a has integer eigenvalues on any irreducible 
representation, and hence from Eq. ITBl 

ol -> 

2—Wk is integer. (14) 

cr 

So on any irreducible representation the eigenvalues of H a are (half)- integer. This fact 
tells us that the magnetic roots b defined by 

B = —b-H (15) 
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Figure 1: The lattice of allowed magnetic charges spanned by the coroots for the gauge group SU(3). 



must lie on a lattice generated by coroots b = -% if the magnetic strength B obeys the Dirac 
condition, Eq. El That is, every point b on the magnetic root lattice is a linear combination 
with integer coefficients of the coroots -%. 

° or 

In Fig. [T]we show the coroot lattice for SU(3). The two simple roots 7^2 span the lattice. 
They define, with H in the fundamental representation, the matrices 

^■H = H 12 = \ diag(l,-l,0) (16) 
j 2 -H = H 23 =\ diag(0,l,-l) (17) 

which indeed have half-integer eigenvalues and are the third components of I-spin and U-spin 
respectively. Often it is more convenient to work with these matrices than with the root vectors 
themselves, as the former are truly simple. 

For SU(iV), the simple roots are given by the generalization of I and U spin. The general 
representative is: 

H kM1 = - diag(0, 0, ...0, 1, -1, 0, ...). (18) 
The first non-zero member is on the k-th diagonal entry, and k ranges from 1 to N, with: 

H N , N+1 = H N>1 = ~ diag(-l, 0, 0, 1). (19) 

The sum of these matrices is zero, and usually the first N — 1 are taken as simple roots. It is 
then clear that we can rephrase the Dirac condition as: 

47T 

B = — {n l2 H 12 + .... + njv-i.jv-H'jv.Jv-i), (20) 
9 

where the n's are integers 1 . 

1 This notation is adapted to our notation for the Wilson loop in later sections. 
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This classifies the possible monopoles for all simple classical Lie-algebras, as hypothesized 
in the second paper of reference j^j. 

For the group SU(2) the consequences of the Dirac condition and this hypothesis are simple. 
We have a doublet with, in units ofAir/g, I3 = ±1/2. Then an iso-triplet with in the same units 
^3 = ±1,0. For a spin J (half)-integer multiplet we have the same. Our matrix gB/4n with 
the spin 1/2 multiplet of magnetic charges gives only on integer spin electric charge multiplets 
an integer. So the magnetic group of SO(3) is SU(2). On the other hand the iso-triplet of 
magnetic charges is compatible with any charge multiplet, half integer or integer, and so the 
magnetic group of SU(2) is SO (3). 

More generally, for the gauge group SU(iV) all possible monopoles are multiplets of a 
magnetic group SU(N)/Z(N). The opposite is also true: the gauge group SXJ(N)/Z(N) admits 
monopoles in multiplets of the magnetic group SU(iV). 

In Fig.[T]the lattice is shown for the gauge group SU(3); for the gauge group SU(3)/Z(3) the 
lattice of monopoles will include the additional sublattice generated by the triplet representa- 
tion. This additional sublattice is obtained in a natural way by introducing the 2 hypercharges 
Yfc, k=l,2. Of course they are not uniquely defined. They generate through exponentiation 
exp(i27rYfc) the center-group elements of Z(3). We may for instance choose a set which is at a 
minimal distance (defined as the trace of the square of the matrix) of the center of the Cartan 
algebra: 

Yt = idiag(2,-l,-l) 

Y 2 = ^diag(l,l,-2). (21) 
In terms of the simple root matrices one finds: 

Y l = ^(2H 12 + H 23 ) 

Y 2 = ^(3H 23 + 2H 12 + H 31 ). (22) 

So following in Fig. the steps along the weight lattice to arrive at Y\ one gets the highest 
weight of the triplet representation. Similarly Y 2 is the highest weight of the anti-triplet 
representation. In Appendix B we formulate this relationship for general SU(iV) and the 
generators of its center-group Z(N). Not only are the N — 1 Y matrices an alternative basis 
for the Cartan algebra. More important for us, they are a measure for the strength of the 
Wilson loops needed to observe the monopoles (see section 14. 2 j) . 

For any general classical Lie group it is the "dual" group [3j built from the dual Lie al- 
gebra [TT] that gives the possible multiplets. The precise dual group with the appropriate 
center-group follows from the same considerations as for the SU(iV) case: the larger the orig- 
inal, electric, gauge group, the more stringent the Dirac condition becomes and a smaller 
magnetic group follows. In mathematical terms it is the connectivity of the group and the 
ensuing Z(N) factors. 

In an earlier paper ^H] precisely these hypothetical monopoles were identified with our 
lumps in 3 dimensions. As we supposed the lumps - now monopoles - to be dilute we can 
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compute their effect on Wilson loops. For SU(iV) groups the choice of adjoint representation 
is unquestionably favoured numerically, as simulated by fc-loops for N = 4, 6 by Teper [23 
and for N = 8 in this paper. 

A comment on the nature of the magnetic group is in order. The monopoles, as bound states 
of magnetic gluons, will transform inside a multiplet under some perhaps very complicated 
function of the original vector potentials. So the global magnetic SU(iV) group will not coincide 
with the original global colour group. This ties in with a phenomenon discovered by the authors 
in Ref. [S|: global colour is not defined on the quantized version of the non-Abelian monopoles, 
due to the long range nature of the colour magnetic fields of the monopole 2 . 



3.1 Monopoles as a dilute gas: the broken symmetry case in the 
Georgi-Glashow model 

At this stage it is useful to put our model into a well-known context, the Georgi-Glashow 
model, with gauge group SU(2) in D = 3 with gauge coupling g 3 . According to our hypothesis 
we have a dilute gas of iso-triplet monopoles which describes the behaviour of Wilson loops. 
Their density is proportional to g®, the only scale before breaking the symmetry. And their 
screening mass is proportional to g\. Adding an adjoint Higgs scalar with a "heavy" VEV v, 
i.e. v 3> #3, will give us the broken phase with heavy 't Hooft-Polyakov monopoles: 

SO (3) -> 17(1). (23) 

This model was studied by semi-classical methods in a seminal paper by Polyakov [Tlj, in 
the limit that g$/v is small. In that limit the diluteness of the monopoles is a fact. The Wilson 
loop tension a is exponentially small, like the density of the monopoles and the screening mass. 
The exponent is on the order of exp — jv/g^, 7 some numerical constant. The result for the 
string tension can be expressed in terms of the density of monopoles um and the magnetic 
screening mass M by combining: 



a 



n M 




and 



into 



16n 



M 



M 




(24) 
(25) 

(26) 



The correction to the tension is due to the authors in Ref. . 

Eqs. and are typical for weak coupling plasmas. The dimensionless ratio of tension 
over screening is proportional to the number of monopoles inside a sphere of radius the screening 
length: 

_ . , r / i r \ ' 

(27) 



a 



16 



M 3 




2 In the work by Bais and collaborators [H] an interesting interpretation of the magnetic group is proposed 
but its discussion falls beyond the scope of this paper. 
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From Eq. |2HJ this number is seen to be large: 



W'^W^Ww^ 1 (28) 

In our model for the strongly interacting symmetric phase (see Eq. |B] and below), this 
very same ratio is small! Physically, what happens is that the strong coupling creates a 
bound state of the original semi-classical monopoles within the magnetic screening radius. 
And indeed, as stated before, Monte-Carlo simulations in the symmetric phase give for the 
ratio ^ 0.046(2) [23j - 



3.2 Broken symmetry: SU(3)/Z(3) and higher groups 

Our next example is the gauge group SU(3)/Z(3) broken by the adjoint to U(l) 2 or to SU(2) x 
U(l)/Z(2). 

The first case is shown in Fig. |21 Every point represents an 't Hooft-Polyakov monopole in 
the corresponding SU(2) subgroup, as in Eq. The Dirac condition carries integers which 
are the topological winding numbers of the Higgs field [Hj. So this case does not go beyond 
what we already knew. As we wil see in the sequel this phase is not realized in simulations, so 
we will not consider this phase anymore. 
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(a) (b) 

Figure 2: The lattice of charges allowed by the quantisation condition, spanned by the simple coroots. In this 
figure also the direction of the Higgs field <f>o in the Cartan subalgebra is indicated. In (a) the stable charges 
for an arbitrary non-degenerate orientation of the Higgs field are indicated by black dots. In that case the 
residual gauge group is ?7(1) x U(l) and all allowed charges correspond to a winding number. In (b) the Higgs 
field is degenerate and leaves the non-Abelian group U (2) unbroken. Now only one component of the magnetic 
charge is the winding number, and in each topological sector only the smallest total charge is conserved. The 
points symmetric with respect to the Higgs field are conjugate through the unbroken group. 



More interesting is the breaking pattern with unbroken group U(2). This model is very 
often used ^Hj for investigations for non-Abelian monopoles. For momenta p ^> v the broken 
phase is perturbative, for momenta much smaller than g\ the coupling becomes strong. We 
expect screening at those distances, including screening of monopoles. 

If we try to construct an 't Hooft-Polyakov monopole in the unbroken SU(2) group along 
7i in Fig. EJ we will fail because the VEV is lacking in that subgroup. Along the root 72 the 
VEV is non-zero, so along that direction the integers still correspond to a winding number. 
Similarly along the direction 71 + 72, obtained by reflection of 72 w.r.t. the direction of the 
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Figure 3: The phase diagram of the 3d SU(3) + adjoint Higgs theory [13]. The open symbols 
are results from the simulations, and the filled circle is the perturbative result [4"4"1 I45j . 

Higgs breaking $o- Their long range magnetic fields are transforming into each other by the 
unbroken gauge transformations. When trying to quantize these monopoles this property poses 
a problem [S] of consistency, which is related to the fact that for the quantized solutions we 
expect screening. The long range field is unstable. 

The mass of these objects is growing with the size of the VEV in the classical approximation. 
In what follows we will assume that this property survives the non-perturbative quantization. 
In terms of our model, the dilute gas of light octet monopoles in the symmetric phase will after 
breaking leave the expected isotriplet in the unbroken SU(2), and two heavy iso-doublets. The 
iso-triplet stays light after breaking, with a density . The iso-doublets are the monopoles we 
described above and live at the lattice points 72 and 71 + 72 in Fig. El Due to their large mass 
they have an exponentially small density like the 't Hooft-Polyakov monopoles in the previous 
example. Note that in both examples the unbroken group defines a neutral singlet, that was 
present before breaking, but has disappeared in the transition between the two phases. 

For illustration we show in Fig. |3] the phase diagram of the electrostatic theory given by 
the action in Eq. |21 for SU(3) by numerical simulation [IS]. The relevant variables are the 
dimensionless combinations x and xy discussed in section 12.11 There is a first order transition 
for small x, that is semi-classically calculable |4"4"1 l4*H] . It marks the border of the region where 
SU(3) symmetry is broken to SU(2) x U(l) and where the global R symmetry is spontaneously 
broken \V\\ . For larger x the transition becomes second order. Above the border there is the 
unbroken phase. This unbroken phase can be smoothly accessed from the broken phase by the 
dotted, second order transition line. It means that putative monopoles in the unbroken phases 
are smooth deformations of the monopoles in the broken phase. 
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Note the absence of a phase where U(l) x U(l) is unbroken but all other generators are 
broken |45j . In general, for SU(iV), the adjoint Higgs does admit for breakings of the type 

su(fc) x su(iv - k) x u(i) mi. 

This phase diagram not only relates the putative monopoles in the unbroken phase to their 
more familiar analogues in the broken phase ^Zj. In order to detect the monopoles one needs 
an operator that measures their flux. This operator is the Stokes version of the spatial Wilson 
loop, and is intimately related to a similar operator for the broken phase. This is the subject 
of the next section. 



4 Flux representation for spatial 't Hooft and Wilson 
loops 

The monopoles have an effect on spatial Wilson loops, because the loops record magnetic flux. 
The traditional representation of the loop as a line integral is not appropriate to quantify the 
effect, and we have to find a flux representation for the loop. For a loop in the U(l) case we 
have Stokes' theorem: 

exp ig j> dl- A = expig JdS-B for U(l). (29) 

Now the non-Abelian case. For a certain class of irreducible representations of SU(iV) 
one finds a simple and useful generalization of the Abelian case. It is due to Diakonov and 
Petrov |1H] • If R is any one of the fully anti-symmetric irreducible representations given by the 
one column Young tableau with k entries it has highest weight Yk (see appendix A; recall the 
property exp (i27iYk) = e t2lTk / N ). Then one finds for the Wilson loop Wr(L): 

TiVexpig <j> dl- A R = J DQ expig J dS ■ Tr {QY k Q^B} . (30) 

The integration is over regular gauge transforms Q. 

In this section the physical ideas behind this Stokes law will be expounded. There are many 
papers jSHj concerning its derivation, but we have not seen any exploring the significance of the 
class of gauge transformations involved, nor the special role played by the fully anti-symmetric 
irreducible representations. First we will make the Stokes theorem plausible by recalling some 
known features |H] of colour electric analogue of the spatial Wilson loop: the spatial 't Hooft 
loop. 



4.1 Flux representation of the spatial 't Hooft loop 

Colour electric flux is confined inside glueballs. It is only above the critical temperature that 
it becomes visible through the area law obeyed by the thermal average of the spatial 't Hooft 
loop. 

The 't Hooft loop is defined as a loop of a Dirac vortex, with strength Zk = e ^/N in the 
center group. The vortex is created by a gauge transformation Q k with a discontinuity z^, 
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when circumnavigating the vortex. The locus of the discontinuity is a surface S spanned by 
the vortex L. 

We take the simplest gauge transformation that does have a discontinuity of this type. 
Gauge transformations are generated by Gauss' operator D ■ E. If 9(S) makes a unit jump 
when going through the surface in the direction of the normal n, then 

V k {L) = expiry" dxTr [Y k {E ■ 3)6 (S)} (31) 

has the required discontinuity. 

On physical states, only the gradient in the covariant derivative D = d + ig[A, counts [4~T] . 
The gluon charge gf a b c A.E is continuous through the surface. So the spatial 't Hooft loop 
becomes: 

Vfc(L) = expi — / dS ■ Tr (on physical states). (32) 

This operator does not look gauge invariant, although on the physical subspace it is. In order 
to bring it in a manifestly gauge-invariant form, we multiply it on the left with a regular gauge 
transformation Q, and on the right with Q' ; a matrix element of between two physical states 
is not affected by this operation. After integration over all regular transformations, 

V k (L) = J DQ expi— J dS ■ Tr {nY k Q^E} (on physical states). (33) 

is manifestly gauge-invariant. 

— # — * 

It seems plausible to obtain its magnetic analogue by replacing E by B and the coupling 
a = (? 2 /47r by a" 1 . That gives the formula for the Wilson loop, Eq. EH 



4.2 The flux representation for the Wilson loop 

The plausibility argument from the preceding section gives a formula which is consistent with 
the expression given in Ref. [49J for Wilson loop in any representation R, with highest weight 
Hr. Let Q be any gauge transformation that is periodic on the loop. Then, with = 
(d-igA)9) and Vfi = dVt + igttA: 

W R (L) = j DVL exp ig J dS ■ Tr j [H R UlBV) - ^-WVt x Vfi^ j . (34) 

This result, proved in Appendix B, differs from that of the plausibility argument through the 
presence of the second term. This term reduces in the SU(2) case to the familiar 't Hooft 
source term. If we limit ourselves to regular gauge transformations, this second term would 
drop out in the equations of motion. 

In the light of this we feel it is justified to make the following assumption: in our model with 
a dilute monopole gas the contribution of the second term is negligible. A second simplification 
occurs when we are only interested in its (thermal) expectation value. This is because in this 
case Wr acts on the left and on the right only on physical states , so the effect of the regular 
gauge transforms Q is undone. 
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There is a further comment related to this Stokes formula . It is derived under the assump- 
tion (see Appendix B) that it is regulated by the SU(iV) asymmetric top The question is 
whether the pure Yang-Mills theory average can be provided with such a regulator. For N = 2 
and in three dimensions the answer is affirmative [H] by adding an adjoint Higgs system and 
letting the VEV go to zero, followed by decoupling the Higgs in the infinite mass limit. The 
VEV is the moment of inertia of the symmetric top. We can not accommodate the extra 
parameter of the asymmetric top, and this is the reason that the Stokes formula is then only 
valid for the fully antisymmetric irreducible representations with highest weight Hr = Y k . The 
reason for this is that the second order Casimir operator takes its minimal value - with fixed 
N-ality k - in the fully antisymmetric representation. 

For general N the answer is analogous, but Nature realizes only a limited set of Higgs 
phases with only one adjoint Higgs field. They are limited to breakings of the type where 
SU(Jfe) x SU(iV - k) x U(l) is still unbroken, k < [N/2] [46 . That implies once more that the 
proof of the Stokes formula is only valid for those highest weights that have this symmetry, 
i.e. of the form Hr = lY k , I a positive integer. For I = 1 this is the weight of the totally 
antisymmetric Young tableaux with k boxes. Appendix B shows that I > 1 is excluded. This 
ends our discussion of Eq. EH 



4.3 Electric flux loop and its expectation value 

Let us return to the electric flux loop, Eq.|S21 In this case the Q integration drops out when one 
acts with \4(L) on a physical state, because the only effect of V k {L) is to multiply intersecting 
Wilson loops in the physical state with a center-group factor (see [H] for more details): 

exp(i27rYfc) = exp(ik2n/N). (35) 

The thermal expectation value of the 't Hooft loop has been calculated analytically at 
high temperature in powers of g{T), including g{T) 3 . This is possible because the effective 
potential is in low orders built up by hard modes (O(T)) and soft modes (O(gT)). The ultra- 
soft magnetic modes come in at higher orders. This potential has a Z(N) symmetry, and the 
thermal expectation value of the loop tension, 

(V k {L))=exp{-p k (T)A(L)), (36) 

is obtained from the tunneling between two vacua, one corresponding to k = and one corre- 
sponding to k. One finds then [20] that: 

Pk (T) = Pi(T)^^, (37) 
up and including two loop order. Concretely, in one loop order the value of p%(T) is [T8] : 

P ^ = A {N - 1)T2 - (38) 

In three loop the above Casimir scaling is slightly invalidated, as it is found to be in lattice 
simulations (23 EH! 
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4.4 Dilute gas approximation for both electric and magnetic flux 
loops 

From the formulae in the preceding subsection one can easily find the behaviour of the tensions 
in terms of k, once one assumes a dilute gas of gluons for the electric loop and a dilute gas 
of monopoles for the magnetic loop at high temperature. A dilute gas of gluons at high 
temperature will certainly disorder the electric loop. The reason is that the flux from one 
single species of gluons is going through the loop only when within screening distance Ie from 
the loop: 

1 / 3 1 , , 

m D y g z JS 1 

Thus all the flux through the loop can only come from the gluon being in a slab of thickness 
Ie and area A(L) of the loop. Thus the flux is approximated by a theta-f unction in the distance 
d from the loop. 

The total flux from a charged gluon is ±1, as follows from the adjoint representation of Y^. 
Thus the height of the theta-function is ±|, because half of the flux is lost on the loop. Its 
effect on the loop is that it picks up a factor 



Vfc(L) = exp 



rl- { ±- 



■1. (40) 



Note that not the value of the charge, but only the multiplicity of the charge depends on k. 
This multiplicity is k(N — k) for each value. 

Now the distribution function of say £ gluons of a given charged species in such a slab 
is peaked around £, the mean number of gluons in the box. Its width should, according to 
thermodynamics, be proportional to £, like e.g. the Poisson distribution: 

= e ^Y- (41) 

The average of the loop is therefore: 

(V k (L))\ one species = = e~ 2i . (42) 

i 

Together with Eq. EE] this means that a single charged gluon species will determine the thermal 
average of the loop to be an area law: 

A{L) Pk \ oae species = 2£ = 2A{L)l E Tl(T). (43) 

Note the absence of k dependence in the outcome! What counts is that the charge is non-zero, 
but its sign is irrelevant 3 . 



3 The reader might be alarmed by our cavalier treatment of the screening of the flux. The flux that 
a gluon at distance d shines through the loop is exponential in d, not a theta-function! One can correct for 
this by dividing the space above and below the loop in parallel slabs of infinitesimal thickness. This means the 
summand in Eq. 021 is replaced by an integral J d(d)exp[l — cos{27r$(d)}]. As a result the factor 2 in Eg. 1431 
increases by a factor 1.64282.... 
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Thus the only way the k dependence comes in is when we take all the charged gluons into 
account. This number, the multiplicity with respect to the charge Y k , is for the adjoint gluon 
multiplet equal to 2k(N — k). It is the number of non-zero entries in the diagonal adjoint 
representation of Y k . We supposed the gluons to be independent; it follows that 

Pfclall species = 2t ■ 2k{N - k) = 41 E n(T) k{N - k) . (44) 

Thus the fc-loop is proportional to the multiplicity of charged gluons with respect to the charge 
Yfe. Note also that equation ESI together with the density of the gluons being ~ T 3 , makes the 
outcome of the one species calculation parametrically identical to the analytic result in Eq. EH1 
The calculation of the magnetic loop is identical. The unit of magnetic charge is Ait/g 
instead of g, but this is canceled by changing from electric to magnetic loop. It is useful to 
realize that the surface integral J dS ■ B for a single magnetic quasiparticle is given by \B, 
where B is the magnetic charge matrix satisfying the Dirac condition (see section (HJ); this 
condition thus directly leads to the phase ir necessary to disorder the Wilson loop, for any 
member of the adjoint representation. So the thermal expectation for the magnetic fc-tension 
is, as for its electric counterpart: 

a k ~l M n M (T) k(N-k). (45) 

Though superficially very alike, there is an important difference between the two tensions in 
units of the respective screening lengths. The electric tension in those units becomes 

l%Pk ~ l%n. (46) 

On the right hand side we have a large number, 0(g~ 3 (T)), for high T. This is the plasma 
condition. It says that an electric screening volume contains a large number of almost free 
gluons. And corrections are in terms of inverse fractional powers of this ratio, as discussed 
already in section 13.11 On the other hand the Wilson k-tension equals: 

l M<*k ~ lli^M- (47) 

Both the magnetic screening and the magnetic density are 0(g 2 T). So in the ratio the coupling 
drops out. Lattice data tell us the l.h.s. is small for all N. The corrections are discussed in 
section 14.61 

For large N the dimensionless quantity timIm is of order 1/N. This is so because the 
magnetic screening length Im can be shown to be given by the ++ mass of the Hamiltonian of 
2+1 dimensional Yang-Mills theory, and therefore behaves parametrically like l/g 2 NT. The 
density of a single monopole species n M should be j^(g 2 NT) 3 , in order to recover a tension of 
0(1). 

4.5 Monopole multiplets other than the adjoint 

Now we use a general multiplet R carrying a unitary representation D^of the magnetic group 
as the magnetic quasi-particles in our model [21]. Its dimension is cLr. The Lie- representative 
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of the charge Y k is written as (Y k ) R and the corresponding group element as D R (Y k ). Of course 
D R {Y k ) = expi{Y k ) R 

As in the previous subsection, the quasi-particle model produces for a given member r 
(r = 1, 2, . . . , d R ) of the multiplet R an area law for the /c-loop with charge Y k , Eq. EH 

W k (L)\ T = exp [-(1 - Re D R (Y k ) TtT )f] . (48) 

The fc-dependence of the tension due to all members of the multiplet is then proportional to: 

o-fc = [d R - ReTr {D R (Y k ) R } = d R - ReTr {exp(iir(Y k ) R )}) l M n M - (49) 

This result is invariant under a gauge rotation, 

Y k -> ttY k rt (50) 

as Eq. EH suggests. And it reduces to 2k (N — k) for the case where R is the adjoint. The reason 
is that (Yfc)adjoint has either or ±1 on the diagonal as argued already in the previous section. 
Hence the formula counts the multiplicity of charged members of the adjoint multiplet. 
For the spinor representation one finds from Eq. 021 the result quoted in Ref. |2Tj : 

a k ~ [N -k cos (((N -k)n/N)) - (N — k) cos (kn/N))] . (51) 

Both adjoint and spinor multiplets are compared to the lattice data in the section on data 
analysis. 



4.6 Corrections 

There are two sources of corrections to the Casimir scaling formula. One is the diluteness, and 
the other are the effects of Bose- Einstein statistics. 

• The diluteness 5 = I^Um = °"i/ m M * s small (~ 0.05, as discussed in section l2~2~J) but 
produces corrections. The use of classical Boltzmann statistics is allowed at large T, since 
the thermal de Broglie wave length 1/T is much smaller than the inter-particle distance 

i/g ( T)T. 

• As we descend in temperature the diluteness stays constant, since we know from the 
results by Laine and Schroeder [HH] that magnetic quantities are determined to a very 
good approximation through all of the plasma phase by their value at very large T in 
3d Yang-Mills theory, and the running of the coupling due to hard radiative corrections. 
Below T c <j\ is virtually constant [HO] • Unfortunately the behaviour of the magnetic mass 
is not known in the cold phase, but we know its value at T = leading to a diluteness 
~ 0.09( see section l2~2"|) . which suggests that it is small at all temperatures. 

• What changes as we go down in T is the ratio of thermal wave length to inter-particle 
distance. So Bose-Einstein statistics kicks in at temperatures on the order of 4T C , where 
g 2 {T) = 0(1). It seems natural that the transition is where Bose-Einstein condensation 
starts. 

In principle the effects due to the small but non-zero diluteness can be computed. Comparison 
to lattice data [213 in the deconfined phase shows that they should be small, on the order of a 
few percent at most. 
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5 On strings in higher representations and ^ corrections 



We now leave the discussion of the adjoint-monopole-gas model and discuss the properties of 
fc-strings from the point of view of the large- N expansion. 

Standard arguments on large N SU(iV) gauge theory [HU E2] based on the planarity of 
Feynman diagrams and the (assumed) confinement of color, imply that gauge invariant states 
have masses of order N° + iV~ 2 , with a width of order iV~ 2 . Also, no bound state of color 
singlet constituents survives the large N limit: the theory is expected to be a theory of free 
'hadrons'. 

It is interesting to consider, at large but finite number of colors, precisely those states 
whose wavefunction contains a significant component which is a direct product of color singlet 
pieces. Phenomenology provides a number of potential examples. A classic example would 
be the deuteron, a very loosely bound state lying only a few MeV under the nucleon-nucleon 
threshold. Another interesting though less firmly established case is the /o(980) meson, whose 
wavefunction has been discussed in terms of a mixture of a kaon-kaon molecule and a four- 
quark state and ref. therein). Again, the state is only a few MeV under the two-kaon 
threshold. 

At a more theoretical level, there are examples in the pure SU(iV) gauge theory in three and 
four dimensions. Consider the theory defined on a finite (but large: L ^> 1/T C ) hypertorus. In 
addition to glueballs, the spectrum contains 'torelon' states (whose mass we denote by rrik(L)) 
which transform non-trivially under the centre symmetry Zjy. The sectors of different A^-ality 
are protected by this global symmetry. Thus for N > A, one may ask whether two fundamental 
torelons can form a bound state lying under the threshold 2nrik = i(L). That this is indeed the 
case was first numerically demonstrated from first principles in the work [23] . Further, at large 
L the states are string-like and one can ask what the ratios of their string tensions are (we 
may use the fundamental l k = V string tension as the reference). An alternative formulation 
of the problem would consider the strings to be open and attached to static sources in the 
appropriate representation [53]. 

For simplicity, we now focus on the k = 2 sector; for N > 4, the screening of the string is 
forbidden by the centre symmetry. In our view, the first question to settle in the context of 
the large- N expansion is, 'What is the 1/N power of the leading correction to the planar limit 
result o"2 = 2<7i?'. Since mi{L) lies under the threshold 2rrik=i(L) at all L j25J, the question 
arises whether one should think of the k = 2 torelon as a weakly bound state of two k — 1 
torelons, or if the colour structure gets completely rearranged into a single 'unfactorisable' color 
singlet piece. In the nucleon-nucleon system, the analogous question is whether the deuteron 
is primarily a bound state of two nucleons, or a 6-quark state. In the first case, considered 
in ISH], the long-distance attractive force between the two k = 1 strings will be driven by the 
exchange of the lightest (0 ++ ) glueball, while the short distance force is essentially given by 
two-gluon exchange. Both effects are indeed [3JJ] suppressed by 1/N 2 with respect to the free 
propagation of two k = 1 strings. Regarding the second configuration, the simplest classical 
string configuration is that of a single string winding twice around a cycle of the hypertorus. 
At large N , the energy of such a configuration is expected to approach threshold from below 
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at a 1/N 2 rate. 

At finite N and asymptotically large L however, we are in presence of two almost degenerate 
configurations lying near threshold. It is therefore imperative to consider the mixing effects 
between these two configurations. To keep the discussion as simple as possible, we may keep the 
transverse spatial dimensions L± finite, so as to separate two-torelon 'scattering states' from 
the weakly bound states we discuss by a finite, fixed gap. As N is increased, this transverse 
volume can be increased as well without affecting the validity of our treatment of the k = 2 
sector as a 2-level system. 

We suppose, following jHlj, that the Hamiltonian of the SU(iV) gauge theory can be ex- 
panded in inverse powers of 1/N: 

H^N) = t^. (52) 

fc=0 

The existence of the t'Hooft limit implies that H a (L) has the same eigenvalues as the Hamil- 
tonian of the SU(oo) theory in the same spatial volume. We consider now the k = 2 flux tubes 
winding around a cycle of the torus as a quantum mechanical two-state system, as was done 
in [33] for the case of the scalar glueball - adjoint Polyakov loop system in intermediate volume. 
Consider on the one hand the state made of two k = 1 non-interacting closed fundamental 
strings, and on the other a single fundamental string with winding number 2; in this basis H Q 
reads H a = 2m k=1 l 2X 2- 

The 'perturbation' describes the deviations from the planar limit. On the diagonal, the 
corrections are 0(1/N 2 ). Indeed, the attractive potential between two fundamental torelons is 
suppressed by the product of two 3- vertices each of which carries a 1/N factor. On the other 
hand, the amplitude of the transition from one of our basis states to the other only contains 
one such vertex, and therefore the off-diagonal element of our 2x2 hamiltonian is 0(1/N). 
The perturbation hamiltonian in our basis reads: 

AH-(~y N2 h/N ) (53) 
^ n ~[ h/N -h 2 /N 2 ) [b6) 

with hi, h 2 and h of order iV°. It is clear that to leading order in 1/N, the resulting energy 
eigenstates are now the symmetric and anti-symmetric linear combinations of our basis states. 
The associated energies are E A = 2m k=1 - £ + 0(1/N 2 ) and E s = 2m k=1 + ± + 0(1/N 2 ). 
We thus reached the perhaps surprising conclusion that the corrections to the mass of the 
lightest k = 2 string are of order 1/N. There is one state below threshold and one above, 
situated symmetrically about the threshold energy, up to 0(1 /N) corrections. We note that h 
(as well as the hi) is expected to grow proportionally to L at large L, since the breaking of the 
string can occur at any point along the string, so that the ratio in the torelon masses directly 
translates into the ratio of the string tensions. 

A caveat particularly relevant to Monte-Carlo simulations is that in all the considerations 
above we have supposed the strings to be long enough to be able to identify the ratios of string 
tensions ratios with the ratio of torelon masses. Since the 1/L string correction ^Tj\ lowers the 
energy of the string, it induces a repulsive force between two fundamental torelons at finite 
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L [32]. Since the binding energy of the strings reduces at large N , L must be increased so 
that the condition 

o x L 2 > N (54) 

is satisfied to ensure that the ratio of torelon loops yields the correct iV-dependence of the 
string tension ratios. If the large N limit is taken at finite L, the ratio of k = 2 to k = 1 
torelon masses will approach 2 with 1/N 2 corrections, because a mixing amplitude only affects 
the energy spectrum at leading order if the 'unperturbed' states are degenerate (to that order). 
And indeed, the two classical-string configurations we took as unperturbed states have different 
1/L corrections: if 71 is the Liischer coefficient of the fundamental string, the direct-product 
configuration of two k = 1 strings has a 2 r y%/L correction, while the fundamental string with 
winding number 2 admits a 71 /2L string correction. As all lattice simulations so far [23123123, 
ours are done in the regime N < a\L 2 < N 2 . The second inequality implies that the energy 
gap (hi — fi2)/N 2 is parametrically smaller than the string vibrational excitations. 

Eq. 1541 however also implies that the vibrational excitations of the strings are separated by 
47r/L gaps which are parametrically much smaller than the mixing energy h/N. We have thus 
neglected the matrix elements of the Hamiltonian Hi between the two states that we focused 
on and the vibrationnally excited states, since we reduced the diagonalisation problem from 
the full Hilbert space to the space spanned by these two states. Although the neglected matrix 
elements can modify the splitting pattern around the threshold, the mixing between the string 
ground states will be enhanced relatively to the other mixings if the condition o\L 2 <C N 2 
holds. In any case, the parametric size of these matrix elements is aL/N and so we must 
expect the corrections to the k = 2 string tension to be of that order. 

In summary, the predictions of the two-state mixing model are: 

1. the energy eigenstates are the anti-symmetric and the symmetric linear combinations of 
the direct-product configuration of two k = 1 strings and the fundamental string with 
winding number 2. 

2. they are split symmetrically around the threshold energy 2m! (up to 0(1/N) corrections). 

3. the splitting energy itself is of order 1/N. 

5.1 Static potentials 

If one considers open strings attached to static 'quarks', the argument takes a slightly different 
form. The relevant quantities here are the static potentials between colour sources in irreducible 
representations of SU(iV). 

The k=l string binds a quark and a distant antiquark together. Similarly the k=2 configu- 
ration can be viewed as two (weakly interacting) strings each joining one of the quarks to one 
of the antiquarks. If we number the quarks by 1 and 2, and the antiquarks by 1 and 2, then 
there are two classical string configurations which are exactly degenerate: the configuration 
where 1 is attached by a string to 1 and 2 is attached to 2, and the other where 1 is attached 
to 2 and 2 to 1. However, the interaction between the strings can take one configuration into 
the other. Therefore a splitting occurs between the symmetric and anti-symmetric linear com- 
binations, corresponding to the static potential splitting between the k = 2 symmetric and 
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anti-symmetric irreducible representations of SU(iV). There is however general agreement that 
screening of the static sources through virtual gluons implies that the string tension obtained 
in either representation at large enough separations is the same; although it can be difficult to 
demonstrate this in Monte-Carlo simulations. 

5.2 A caveat on the implications of factorisation at large N 

The standard way to extract the static potential for fundamental charges, namely by measuring 
the expectation value of a rectangular Wilson loop of size RxT, T 3> R, can be generalised to 
extract it for any representation jnSj. In particular, the simplest way to obtain a representation 
of 7V-ality k = 2 is to take the real part of the square of W(R, T), the trace of the fundamental 
Wilson loop. At finite R, T, the factorisation property of gauge invariant operators (see for 
instance [56J) then implies that the expectation value of this operator is given by 

(W(R,T) 2 } = (W(R,T)} 2 (1 + 0{l/N 2 )). (factorisation) (55) 

On the other hand, if we consider small separations R, asymptotic freedom implies that the 
short-distance potential in an irreducible representation 1Z is given by C-R,a s /R. The symmetric 
k = 2 representation has C-r — Cs — 2(N + 2)Cp/(N + 1), while the k = 2 anti-symmetric has 
Ca = 2(iV — 2)Cf/{N — 1). In particular, for the fundamental representation, it is 

W(R,T) = exp (-^j + 0(1/N 2 ), R « a^ 2 . (56) 

The operator W 2 (R, T) belongs to a representation that can be reduced into the symmetric 
and anti-symmetric. Therefore, if we take the T — > oo limit, the potential energy of the 
anti-symmetric representation dominates the expectation value of W 2 (R, T): 

lim {W 2 (R, T)> oc e -CA«.CT/R = T ^2 e f£ (1 + (1/N 2 )), R < a' 1/2 . (57) 

Thus tree-level perturbation theory contradicts the large- N counting rules concerning the lead- 
ing corrections to factorisation. The origin of the paradox lies in the straightforward T — > oo 
limit necessary to filter out the ground state. If the contribution from the symmetric repre- 
sentation is kept, the large N limit of the small- R, large-T Wilson loop (W 2 (R, T)) is given 
by 

hm (W 2 (R,T)) cx (W(R, T)) 2 cosh (1 + 0(1/N 2 )) } (58) 

which, at fixed T, has 1/N 2 corrections to the planar limit result. 

What have we learnt? The large- N factorisation property does not necessarily imply that 
the lowest energy state of a 'meson' made of a static colour source in a certain representation 
and its anti-source has 0(1/ N 2 ) corrections, because other representations of same M-ality 
become degenerate with it in the N —>■ oo limit. Since string tensions are extracted from the 
lowest energy at large R, the same caveat applies to them. 
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5.3 Strings in open and closed form 

We finish with a remark on the relation between different representations of same 7V-ality and 
excited states in the open and the closed string sectors. To that end it is useful to consider the 
correlator of Polyakov loops of length L (P-^(0)P-^/(x)}. This expectation value is interpreted 
(from the point of view of the transfer matrix along the dimension of size L) as the free 
energy of the system in the presence of two static charges in the given representations. When 
1Z — TV — N, the fundamental representation, the heavy-heavy bound state can a priori be in 
the adjoint or the singlet representation (in SU(3): 3 ® 3 = 8 © 1). Now, it is believed that 
only bound states in the singlet representation have a finite free-energy in the confined phase. 
That means that if the heavy charges themselves are not in the singlet representation, virtual 
gluons will try and screen the chromoelectric field emanating from this coloured bound state 
until it is a singlet again. Since the gluons are in the adjoint representation, they can screen 
the configurations of the heavy-quark bound state that are in the adjoint representation, albeit 
at a certain energy cost. On the other hand they cannot screen a single heavy quark, and the 
latter therefore has an infinite free energy. 

Suppose we want to determine the static potential for sources in all possible representations 
of SU(iV) (not necessarily irreducible) of a given A/"-ality k and up to a given size. Clearly it 
is sufficient to determine the Polyakov loop correlators between the irreducible representations 
obtained in the decomposition of the direct product representation of k quarks. The question 
then arises whether the cross-correlations (i.e. for 1Z ^ 1Z') between Polyakov loops in these 
irreducible representations vanish or not. If they do, it implies that the energy eigenstates are 
in definite irreducible representations. 

Consider the k = 2 case. The direct product of two fundamental representations decompose 
into a symmetric and an anti-symmetric representation: iV ® N = A © 5*. In the most familiar 
case of SU(3), the anti-symmetric representation is nothing but the 3 (anti-fundamental): 
3 © 3 = 6 + 3. So we are asking whether (Ps(0)P^(x)) has to vanish. We have 6 © 3 = 10 © 8, 
so that virtual gluons can screen the adjoint piece, thus ensuring that the free energy of 
this system is finite. So in general these cross- correlations do not vanish. It is easy to see 
(using Young tableaux) that in SU(iV) the adjoint representation appears exactly once in the 
decomposition of S © A. However since gluons have to screen the heavy-heavy system, the 
(Ps(0)Pa(x)) are 1/N 2 suppressed at large N. Let us now see what this conclusion implies for 
the determination of the fc-string tensions in the open and the closed string sectors. 

In order to study the lightest open string, one may in principle choose to immerse any one 
static source of the relevant 7V-ality in the system, since for large enough R the sources are 
expected to be screened down to the representation with the smallest string tension. Once 
the linear behaviour of V(R) with the latter slope sets in, the differences between the static 
potentials in irreducible representations of same A/"-ality are expected to become only weakly 
independent (they correspond to 'gluelump' masses [Ml)- For long enough strings, the lowest 
excitations of any of these static 'mesons' correspond to the lowest excitations of that string, 
which come in gaps of order 1/R. In short, there is at most one stable open string for a given 
7V-ality. 

It is also possible to interpret the Polyakov loop correlator with a transfer matrix along the 
direction x. One is then measuring the spectrum of states of the gauge theory which carry a 
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winding number with respect to a cycle of the hypertorus of length L. Of course, since the 
Polyakov loop correlator has a unique asymptotic area law, the coefficient in front of the area 
defines both the string tension in the open as in the closed string sector. Just as in the open 
string case, there cannot be more than one stable string per TV-ality because of the screening 
by gluons. A simple picture |2Sj is that virtual gluons screen the unstable string down to the 
stable one and propagate along it until they annihilate around the cycle of the torus. 

For long enough torelons, the lowest closed-string excitations are again expected to be 
string-like, i.e. coming in 1/L gaps. There can be resonant states of the torelons (lying above 
the fc-torelon threshold) whose energies grow linearly with L. It is then natural to associate 
them with meta-stable strings. 

What we inferred about the cross-correlations between different irreducible representations 
above tells us that the energy eigenstates do not in general belong to irreducible representations 
of SU(iV), although the mixing between them is suppressed (at least in the k = 2 case) by 
1/N 2 . 

6 Lattice simulations 

We extract string tensions in the three-dimensional SU(8) gauge theory from the masses of 
'torelons', gauge invariant states transforming non-trivially under the Z(N) symmetry of the 
action; they wind around one spatial cycle of a the hypertorus. These masses are extracted 
from the exponential decay of correlation functions at 'large' Euclidean time. To enhance 
the signal-to-noise ratio, we use fuzzing techniques in the construction of our operators as 
described in [29J. The correlation functions are measured on gauge configurations generated 
by a Monte-Carlo program. We use the original Wilson action jSH]- The configuration is 
updated by sequences of 'sweeps'. One sweep consists of updating all links by performing 
either a heat-bath (HB) jEZ| or an over-relaxation (OR) [63] step on N(N - l)/2 of its SU(2) 
subgroups [HI]. The ratio of HB:OR is 1:3, and we typically perform a sequence of 1 HB and 
3 OR between measurements. We use a 2-level algorithm [32] as described in [33]. The latter 
reference also contains a detailed comparison of efficiency of the ordinary 1-level and 2-level 
algorithms. The number of measurements performed at fixed time-slices was 800 at f3 — 115, 
200 at = 138 and 40 at (3 = 172.5. 

6.1 String corrections 

Consider the Euclidean gauge theory on a L x L x T hypertorus, with cycles of length L. The 
gauge-invariant states with winding number k ^ around one spatial cycle of the hypertorus 
are called torelons. In the Hamiltonian language they are created by spatial Polyakov-loop 
operators with A/"-ality k; a description of the operators used can be found in appendix C. If 
the dynamics of a torelon state of length Ly/a 3> 1 is described by an effective string action, 
then the expression for its mass as a function of its length reads 




(59) 



24 



where 7 is a numerical coefficient of order one which only depends on the universality class 
of the string [T7j. Recent accurate numerical results jS7| 122] show that the flux-tube in the 
fundamental representation belongs to the bosonic string class. In the case of a torelon this 
implies that 

7 = 76 = V - ■ 60 


In general, if Eq. holds, the ratio of the lightest fc-torelon mass to the k = 1 torelon 
mass is given by 

m k(T , a k a k ( 1 \ 2 7i<j fc 

— \L ) = — + + ° > "fc = 7fc 61 

mi cti eriZ/ \GiL 2 J ox 

The sign of a& is of interest. If the fc-string is a weakly bound state of k fundamental strings, 
then one would expect 7^ = kji and therefore a k = — (k — ^)7i < 0. If one the other hand the 
fluctuations of the k fundamental strings are 'in phase', then the number of degrees of freedom 
on the worldsheet of the fc-string is the same as for the fundamental string, hence 7^ = 71 and 
a k = (a - i) 7l > 0. 

At each lattice spacing, we measured the masses of the torelon states of at least three 
different lengths. In practice, we use two asymmetric lattices of the type L\ x L2 x L t and 
L2 x L3 x L t . In this way, we obtain three different lengths of the torelon, and we can also 
check for any dependence on the transverse size of the lattice by comparing the mass obtained 
for the torelon of length L 2 on the two lattices. The Lj range from 1.4fm to 3fm, if we set the 
scale by ^fo\ = 440MeV. This is longer than what has been normally measured so far, and is 
made possible by the use of the two- level algorithm. We use Eq. EH to obtain the fundamental 
string tension by fitting m(L)/L with a linear function in 1/L 2 . The intercept yields the string 
tension; the slope gives the Liischer coefficient. Whether the functional form (JHH|) successfully 
describes the leading deviation from constant linear mass density is controlled by the \ 2 of the 
fit. 

Systematic errors play an important role in comparing the numerical data to model pre- 
dictions. In an attempt to get them under control we propose two separate ways to extract 
the ratios of string tensions (we refer to the first method as the 'unconstrained' one, and the 
second as the 'constrained' one). In practice, having learnt from the pros and contras of both 
data analyses, we present our final, 'educated' analysis in section IB. 2 .41 



1. Firstly the ratios of torelon masses m k (L)/mi(L) are fitted according to Eq. EH with a 
linear function in 1/L 2 , and the intercept gives us the ratio In this way, we need make no 
assumption about the values of the coefficients 7 corresponding to the different representations; 
in particular, the different strings could have different coefficients 7^. Finally, these string ratios 
are extrapolated to the continuum, a — > 0, in a standard way. 



2. The second analysis will assume that all fc-strings belong to the bosonic class. Conse- 
quently, we can extract the string tension ratio from Eq. lousing the estimate a k ~ — 7^. 
with 7^. = 71 = 7fe at every L. The estimates of the ratios obtained at different L are then 
simply averaged, as long as they are compatible with eachother, to produce the estimates of 
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the string-tension ratios. If the \ 2 °f t ne average is large, we drop the smallest L until an 
acceptable \ 2 is reached. The continuum limit is then taken. 

The multi-level algorithm allows us to apply the variational method [HOI on the correlation 
matrices at t > 2a of Euclidean time separation, an improvement over the traditional where 
the method is usually unstable unless t = 0, although the method really finds its justification 
when applied at large t. 

Having said that, we note that this work constitutes the first attempt to extract the k = 4 
string tension from Monte-Carlo simulations, and should be regarded as exploratory in that 
sector. Indeed we found that the variational method jHOj generally became unstable if all 
five operators listed in appendix C were fed in the generalised eigenvalue problem. As a 
consequence only three or four of the five types of measured operators (at the 'best' level of 
smearing-blocking) were finally employed. This and the fact that we only have a short range 
in Euclidean time to identify the mass plateau, due to the rapid fall-off of the signal, means 
that the k = 4 string tension has a significant systematic error attached to it. For the lower 
k states, these problems are less accute and we are much more confident about their mass 
estimates. 

6.2 Data analysis 

We give the masses of the lightest spatial torelons of each A/"-ality in Tab. Q Tab. 0] gives 
estimates of the first-excited torelon mass in the k = 2 sector, that will be discussed below. 
Within the range considered (L ~ 1.9fm, 0.8L < L± < 1.2L), we certainly find no dependence 
of the k = 1 torelon masses on the transverse size. There is also no statistically significant 
variation of the lightest higher-/c torelon masses. Transverse size corrections are expected to 
be suppressed by a power oil/ L varying continuously with L±, but greater than 3 [55| . 

We show on Fig. El the local effective mass of the correlators in the k — 1 and k = 2 
representations. We emphasize that the variational method, which yields (quasi-) orthogonal 
states, automatically picks out the symmetric and anti-symmetric linear combinations (within 
very small fluctuations on the coefficients). We shall come back to this point in the discussion 
below, section I6~3l 

6.2.1 Setting the scale 

Although one could choose the (dimensionful) coupling to set the scale, we prefer to use y/o\ 
for this purpose. We extract the fundamental string tension in lattice units at each of our 
three lattice spacings by linearly extrapolating the torelon mass per unit length, mr/g^L as a 
function of l/(g 4 L 2 ), to infinite L. This is illustrated by Fig. H]in the case (3 = |^ = 138. The 
resulting string tensions are given in Tab. EJ We are able to extract the coefficient of the 1/R 
string correction with moderate accuracy; it is also given in Tab. EJ The coefficients we obtain 
are within 1.3 standard deviations of the bosonic string value. 

Similarly, we can extract the k = 2 string tension and its string correction coefficient 72 
(Fig. bottom plot). It is clear however that the accuracy of the data does not allow us to 
estimate 72 . 
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6.2.2 Unconstrained extrapolations 

In this analysis, for each lattice spacing we extrapolate the ratios of fc-torelon masses to L = oo, 
assuming l/(g 2 L) 2 corrections. In most cases, we have three torelon lengths to extrapolate. 
For the intermediate length, where we have two statistically independent and compatible values 
obtained at different transverse sizes of the spatial lattice, the average (weighted by the inverse 
square of the statistical error) of the two values was taken, whilst keeping the smaller of the two 
errors. In the ratio of the fc-torelon to the k — 1 torelon mass obtained in the same simulation, 
we checked in several cases that the error bars obtained by assuming statistical independence 
do not differ by more than 10% from the jacknife values of the error bars; the former are then 
used in the following. 

We note that the results of these extrapolations done at different lattice spacings are in fact 
consistent within error bars (see Tab. El); it appears that finite lattice spacing effects are much 
smaller than the finite string- length effects in our data set. The \ 2 °f each of these fits are 
good (smaller than 1), except for the extrapolation of the 02/01 at (3 — 172.5, where x 2 — 3-0. 
Since the L = 00 extrapolated value is entirely consistent with that obtained at the other 
values of /3, we attribute this to a statistical fluctuation and, perhaps, a slight underestimation 
of the error bars (due to the neglect of the sort of systematic errors mentioned at the end of 
section 16. lj) . 

Now extrapolating these string tension ratios to the continuum (assuming 0(01 a 2 ) discreti- 
sation errors), we obtain 02/01 = 1.701(77), 03/01 = 2.31(16) and 04/01 = 1.96(23). The \ 2 
of these fits are smaller than 1. The final error bars have blown up due to a somewhat small 
level-arm in the continuum extrapolation. 

6.2.3 Constrained analysis 

In this independent analysis, we assume the validity of Eq.EHlwith 7 given by the bosonic string 
value Eq. inHto extract the string tensions at finite L (neglecting the 0(1/ X 4 ) terms); see the 
string tension ratios in Tab. |H1 where again statistical errors have been added in quadrature. 
In most cases, these ratios are consistent with being independent of L for L > 1.4fm. The 
exceptions concern the k = 2 string at the two smaller lattice spacings (due to the accuracy 
of the data), where we drop the smallest L in our average. We note that, compared to the 
values of the unconstrained analysis (Tab. |2j, the ratios are systematically larger. The ratios 
for the fc-strings in the continuum limit now are: 02/01 = 1.776(33), 03/01 = 2.210(50) and 
04/01 = 2.282(63). The x 2 of the fits are again smaller than 1. 

6.2.4 Final 'educated' analysis 

We consider the preceding analysis to be somewhat unsatifactory, because it assumes a specific 
correction to the /c-string energies which we are not presently able to confirm directly (see 
Fig. EJ, and yet (in the k = 2 case) is of the same order of magnitude as the difference 
between two theoretical expectations we are to compare our data to. Moreover we saw that 
the string tension ratios obtained in this way are systematically higher than if we do not make 
any assumptions about the Luscher coefficients, although the trend is at the one-standard- 
deviation level. 
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The first analysis is well-principled but suffers from the succession of extrapolations to 
L = oo and a = 0, most of which are based on three data points only and are therefore rather 
unstable. Considering the large- L extrapolation (Fig. in particular the k = 2 plot), we see 
that while the coarsest lattice spacing data still shows a difference with respect to the other 
two data sets, the latter two essentially fall on a single curve. Therefore we drop the j3 — 115 
data and combine the data at (3 — 138 and (3 = 172.5 to do a single extrapolation to L = oo. 
The result is: 

lattice (final) adj.monop. fund.monop. trigonometric 
0-2/0-1 = 1.707(28) 1.714 2.105 1.848 

0-3/0! = 2.182(55) 2.143 2.958 2.414 (62) 

04/01 = 2.203(82) 2.286 3.256 2.613. 

(the x 2 are respectively 3.5/5, 2.35/4 and 2.6/4). One ought to associate a systematic error 
with this final result which is of the same order as the statistical error, since evidence for the 
absence of scaling violations was given only at that level of accuracy. We also note that the 
slope, which corresponds to the quantity a& defined in section EJ is clearly positive, clearly 
demonstrating that the central charge of a fc-string is not k times that of the fundamental 
string. 

It is hoped that presenting different analysis strategies has given the reader a sense of the 
challenge presented by these calculations to reduce the systematic errors on the final string 
tension ratios. Comparing our data to the theoretical predictions of various models (Eq. Ifi2*|) 
we find that our data is consistent with the Casimir scaling predicted by the adjoint monopole 
model, and rules out the sine formula by at least 3 standard deviations at all k (even if 
we conservatively assign to the data a systematic error equal to the statistical one). These 
conclusions agree with earlier results obtained for SU(4) and SU(6) [22], although the accuracy 
was high enough for N = 4 to see a (non unexpected) small deviation from Casimir scaling. 



6.3 Excited k = 2 strings 

Fig. El shows the local effective masses, defined as m eS (t + |) = log ( ct+a) ) > °^ severa l °f our 
operators; a plateau is the signature that an energy eigenstate is saturating the correlator. We 
show the local effective mass for our best k — 1 operator. The latter has been determined 
by a variational method [HD] allowing to minimise the contributions from excited states to 
the correlator. Although several levels of fuzzing were included in the variational basis, the 
output wave function turned out to be dominated by a single level of fuzzing. We note that 
its plateau extends out to t ~ -1 ^ 2 , giving as confidence in our mass extraction. In the 
k = 2 sector, we show local effective masses corresponding to the same level of fuzzing that 
was optimal for the k = 1 sector (the inclusion of other fuzzing levels leads to imperceptible 
changes in the mass plateaux). After the basis operators had been normalised in such a way 
that (Oi(0)O*(t = 0)) = 1 (i = 1,2), the variational procedure selected (within 0(1%) error 
bars) the anti-symmetric and the symmetric linear combinations of the operators 0\ = Tr {P 2 } 
and O2 = (Tr P) 2 for respectively the lightest state and the first excited state. Correspondingly 
these operators show quite convincing mass plateaux. By comparison, the individual operators 
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have a less good overlap onto the lightest state, although the signal extends far enough in 
Euclidean time to see that this overlap is not strongly suppressed: their local-effective- masses 
end up being consistent with the plateau of the anti-symmetric combination. Remarkably their 
whole correlators seem to agree at all t. 

Thus the theoretical expectation that the energy eigenstates belong to irreducible repre- 
sentations of SU(iV) up to 0(1/N 2 ) admixtures, which was motivated both in the two-state 
mixing model and by more general arguments about the iV-dependence of screening (resp. 
sections 151 and I5~^|) . is indeed well verified. 

Another prediction of the two-state mixing model presented in section El is that the lightest 
and the first-excited states should be split symmetrically around the threshold energy of 2mk=i 
(to leading order in 1/N). This is tested quantitatively in Tab. EJ which directly compares 
2m k= i to \{m k= 2 + m* k=2 ). The latter two quantities are remarkably close for all string lengths 
and lattice spacings, and in many cases they are compatible within the quite small error bars 4 . 

As we discuss next, the numerical evidence obtained so far favours a binding energy of 
fc-strings of order 1/N . The three predictions that follow straightforwardly from the two-state 
mixing model presented in section have thus been verified quantitatively. 

6.4 iV-dependence of the binding energy of /c-strings 

On Fig. (top) we show the relative binding energy of fc-strings per unit length, ka% — a k , in 
units of <j\ and rescaled by a factor N. We do so by compiling our SU(8) lattice data with the 
SU(4) and SU(6) data from [23]. The predictions of Casimir scaling and of the Sine formula 
are also plotted. The figure certainly suggests that the k = 2 binding energy scales as 1/N, 
with a coefficient of order one. By contrast, to account for the measured N = 8, k = 2 binding 
energy in a 1/N 2 expansion, the first coefficient would have to be about 20. We further note 
that the numerical agreement between the Casimir scaling prediction and the lattice data is 
quite remarkable. If anything, it lies somewhat above the lattice data, indicating that the 
fc-strings are slightly less tightly bound that the Casimir formula suggests. 

On the bottom plot, we show the k = N/2 string tension, rescaled by a factor 2/N, as a 
function of 1/N. The data is plausibly heading towards a finite value at iV = oo. Here too, 
Casimir scaling offers a good description of the data. Note that it predicts that the binding 
energy of the k = N/2 string is half of the energy of k non-interacting fundamental strings. The 
case k = N/2 is special in that the relevant operators (listed in appendix C) and their complex 
conjugate can mix through the appearance of the baryonic vertex on the string. Pictorially it 
swaps the oriented string from on orientation to the other. Naturally the eigenstates of the 
Hamiltonian are also eigenstates of the charge conjugation operator, i.e. the real and imaginary 
parts of the operators, which are respectively C = + and C — —. However the existence of 
a non-vanishing transition probability between the strings of definite orientations means that 
there is a splitting between the C = + and the C = — states of A/"-ality N/2 (for k < N/2, the 
center symmetry forces the degeneracy of these two sets of states). In a two-state Hamiltonian 
formalism, the Casimir formula thus suggests that the Hamiltonian matrix element (per unit 

4 As a technical aside, we note that it is essential here to use the correlations between the local effective mass 
of the lightest and the first-excited k = 2 states, as they seem to be strongly anti-correlated. 
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length of the string) associated with the baryon-vertex is yd. 

It should be noted that the cost of computing the binding energy for a given k naively 
increases as iV 5 (N 3 for the cost of multiplying SU(iV) together in the Monte-Carlo simulation 
and a oc N 2 increase of the statistics to compensate for the 1/N size of the binding energy). 
And this does not even take into account the condition o\I? ^> N formulated in section EJ 
Therefore it could be useful to also compute the string tension ratios for SU(5) and SU(7) 
before moving to even larger groups. 

7 Conclusion 

The picture of hot multi-color QCD considered in this paper relies on the separation of the 
hard, soft and ultra-soft scales by means of a well-known sequence of two effective theories. 
The second of these, which describes the magnetic properties of the quark-gluon plasma at 
energies of order g 2 (T)T, is the 3D SU(iV) gauge theory obtained by dimensionally reducing 
the original theory. 

In the case of the spatial t'Hooft fc-loop, which records the fluctuations of the electric flux 
of A/"-ality k going through it, perturbation theory is directly applicable and becomes ever more 
accurate at higher temperatures. On the other hand, a simple physical picture emerges if the 
electric-flux fluctuations are attributed to gluons passing randomly through it. By assuming 
these quasi-particles to be dilute and non-interacting, one easily derives an expression whose 
parametric dependence on the temperature and the A/"-ality of the loop match the perturbative 
result. 

The adjoint monopole gas model discussed in this paper assumes a similar picture to hold for 
the spatial Wilson /c-loops: the magnetic-flux is attributed to non-Abelian monopoles in the 3D 
SU(iV) theory. The assumption that these monopoles are in the adjoint representation directly 
leads to the prediction of Casimir scaling for the ratios of the associated /c-string tensions. 
The fc-dependence is indeed given by the multiplicity of charged monopoles with respect to 
the charge Yj, of the loop, while the sensitivity to other details of the model is reduced in 
these ratios. Perturbation theory is not applicable in this sector, but non-perturbative lattice 
Monte-Carlo calculations of the fc-string tension ratios, although numerically challenging, are 
in principle straightforward. 

Previous simulations for N = 4, 6, as well as the N = 8 data presented in this paper 
confirm the Casimir scaling property of A;-string tension at the few percent level. Care must 
be taken in these calculations that the strings are long enough, oL 2 ^> N, for the quantum 
corrections to the string energy to be subleading with respect to the weak binding energy of the 
k fundamental strings. The energy of the string is then large and a multi-level algorithm jHZ| 
proved useful in this situation to reduce the variance on the correlator from which this energy 
is extracted. While the Casimir scaling prediction lies slightly above the 3-loop expression 
for the t'Hooft loop, it is slightly lower than the lattice results for spatial Wilson loops: the 
(small) corrections to Casimir scaling in the magnetic sector seem to have the opposite sign 
with respect to the electric sector. 

The fundamental assumptions of the model can be further tested. The behaviour of k- 
loops, as we argued, is the simplest observable to consider. It can be measured for other gauge 
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groups, as long as the center is Z(N), N > 4; indeed the adjoint monopoles may be formed 
in any non-Abelian gauge theory. The classical groups have as candidates, apart from the 
SU(iV) groups, the Spin(4p + 2) groups with center group Z(4). One may also introduce an 
adjoint Higgs field which acquires a VEV. Depending on the symmetry breaking pattern, some 
monopoles will become heavy, and the fc-ratios will change in a predictable way |47j . 

We also discussed the behaviour of the 't Hooft loops for different 7V-alities at finite tem- 
perature. Here the same Casimir scaling is observed in the lattice data [211 12H] as predicted 
by perturbation theory for high T (and by the quasi-particle picture of gluons). Surprisingly, 
the scaling continues to hold down to practically T c . And the same is true for the magnetic 
k-loops [HBIPOJ: Casimir scaling stays valid down to ~ T c . 

In conclusion: all available data on electric and magnetic k-loops for N > 4 are consistent 
with a screened electric and magnetic quasi-particle model throughout all of the plasma phase. 
Whether the same is true for QCD (N = 3) remains to be tested. 
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Appendix A 

In this appendix we briefly indicate the group theory needed to get from a given Young tableau 
(defining the irreducible representation R) the corresponding highest weight and the value of 
the quadratic Casimir. In what follows we suppose a representation to be irreducible without 
mentioning so. We choose the Yk as follows: 

F fc = i- diag( fc, k, . . . , k , k — N, k — N, . . . , k — N ) . (63) 

N—k times k times 
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Let the Young tableau have n\ boxes in the first row, n 2 in the second row, etc. . . Then one 
can define the non-negative numbers wi — n\ — rty + i. Now the highest weight matrix for the 
Young tableau is defined through the Y k matrices: 

JV-l 

H R =Y. w i Y i- (64) 
1=1 

For example, for the totally antisymmetric tableau of k boxes in one column we have Hr = Y k . 
For the totally symmetric tableau with all k boxes in one row Hr = kL\. Note that the stability 
group ofYfc (the subgroup of SU(iV) matrices commuting with Y k ) is SU(fc) x SU(iV — k) x U(l). 
So the totally antisymmetric representation with k squares has a highest weight with this 
stability group. All other representations with k squares have different stability groups. 
We define one more diagonal N x N matrix by: 

N-l 

2Y = 2j2Yi = diag(iV - 1, N - 3, -N + 1). (65) 
i=i 

The quadratic Casimir operator C 2 (N, k, {wi}) = C 2 (R) is defined by summing the square of 
all generators T a in the representation R. The result is YJT^ = C 2 (R)1r, where 1r is the unit 
matrix in R and C 2 {R) is a c-number (normalization is [T a ,Tb] = if a bcT c , fabcfbcd = N5 a d)- 
Then the quadratic Casimir equals: 

C 2 (R) = ±(Tr{H R } + 2Tr{Y H R }). (66) 

The quadratic Casimir for the fundamental representation is CV = ■ The Casimir for 

the antisymmetric representation is then 

C 2 (R=AS) = C F ^^, (67) 



and for the symmetric representation it is: 

C 2 (R = SS)=C F k -^±^. (68) 
To derive these relations one needs the inner product of two Y matrices: 

Tr {Y k Y t } = ^(mm(k, l)N - kl). (69) 

One can show that for fixed k < N the antisymmetric Casimir is the minimal one. 

Finally we give the relation between the Y k matrices and the Chevalley basis H k ,k+i = 
diag(0, 0,0, ....1, —1,0, ...0), Y N>N+ i = Y N>1 . We introduce the matrices y k , with y 1 = NY 1 . y 2 
follows from yi by a cyclic permutation of the diagonal elements: the first becomes the second 
and so forth. y 3 follows from y 2 the same way. We keep doing this until we have reached y^, 
yN+i = yi- The sum of all the y k vanishes. The Y k are related to the y k by: 

NY k = j^ yi . (70) 
i=i 
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(N - 1)H 12 + (N- 2)H 23 + {N- 3)H M + + l.H N . 1N + O.H N1 

0. H 12 +(N- 1)H 23 + (N — 2)H 3A + + 2.HN - IN + l.H m 

1. H 12 + O.H 23 + (N — 1)H U + + 3.HN - IN + 2.H m 

The general term is: 

X -y k — (k- 2)H 12 + (k- 3)H 23 + ... + O.H kk ^ + (N - l)H kk + .... + kH N . 1N + (k — l)H Nl . (71) 

Note the diagonal translation invariance of the coefficient matrix M relating the N y k to the 
N H ljl+1 . 

Evaluate the first diagonal element of the r.h.s. in Eq. EJ using the matrix M above. It 
equals, because of the relative sign in the non-zero elements of H 12 and H m , the difference of 
the first column and the last column of M, up and including the first k rows of M. Because of 
the translation invariance only the difference Mn — M k ^ = N — k survives. The second diagonal 
element equals M k \ — M\ 2 = N — k. This goes on till we reach the coefficient M kk _\ = 0. Then 
there is a jump to M kjk+ i = N — 1, and the diagonal elements become —k. So we reproduced 
the matrix NY k in Eq. E3 The charges NY k lie on the root lattice spanned by the Hu + i. 



Then: 



1 




1 

nV3 = 



Appendix B 

The derivation of Eq. |^] is based on simple properties of the quantum-mechanical SU(iV) 
rotator in an external field. We reproduce it here in a form that should render its origins clear. 

We are interested in the Wilson line between two points x(sx) and x(s 2 ), the line between 
the two points being parametrized by s. The line is the ordered product (from right to left) in 
some irreducible representation R Q with highest weight H° of unitary matrices with dimension 
d Ro : 

(72) 



f s f 

W(sf,Si) = P exp ig I A s ds 

J Si 



Here A s = j^-A the projection of A on the line. 

The Wilson loop is covariantly constant along the curve L, d s W\s, s^) — igW^A s = 0, so 
along the loop one can write the vector potential as a pure gauge 

A s = —Ud s rf (73) 
W 

with U = W(s,Si)Ui with Ui an arbitrary SU(iV) matrix. So the loop becomes in the repre- 
sentation Rq the unitary matrix : 

D Ro (U f Uj). (74) 

We close the Wilson loop (so Sf and Sj represent the same point) and take the trace of Eq. ITU 
The main result of this appendix is: 
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• The sum over all irreducible representations of the normalized trace of the Wilson loop 
is the group average of the propagator of the SU(iV) rotator with Hamiltonian Ti in the 
external field H°, Eq. E21 

• Then a special limit of the path integral version of this propagator, Eq. EH1 is nothing 
but Eq. El 

To define the Hamiltonian we start with the generators l a of the group SU(iV). They act 
as left multiplication on the Hilbert space defined on the group manifold, with kets \U), [/an 
SU(iV) element. The generators obey the commutation relations 

[LJb] = ifabJc, With fabcfa'bc = N8 aa >. (75) 

The N — 1 diagonal generators Id are denoted by the suffix d, whereas a general generator 
is denoted by indices a. 

Likewise we can define the generators r a of right multiplication. 

The state space of kets |0 > is has an complete orthogonal basis consisting of all the 
irreducible representations D R (Q)^ r , and the orthogonality relations read: 

/ dn l>* (Q)< n (fit) = _L 6RRl6kn6l ^ (76) 

The dimension of the representation R equals Tr {D R (1)} = d^. We define kets and bras 
\R; I, r) with the property: 

(n\R;l,r) =D R (Q) l>r . (77) 

From Eq. [23 one sees that the norm of such kets is d^ 1 . 

The Hamiltonian is defined in terms of the left multiplication generators as: 

^ = ^(E^-E^ + E(^-^ ) 2 ). (78) 

The last term represents the coupling of the rotator to the highest weight H° written in 
component form H = J2dHd^, which gives a Zeeman effect for the energy levels in any 
representation R. 

Write J2 a \\ = C*2(-R) for the value of the quadratic Casimir in R and denote I for the N — 1 
diagonal quantum numbers Id in D(Q)i )r . Similarly the N — 1 diagonal right multiplication 
generators have quantum numbers, denoted by r. Then one has for the eigenvalues E° of this 
Hamiltonian: 

e°(R, i) = ^(c 2 (R) -£^) + 27 EC- - \ H ") 2 - ( 79 ) 

d d 

Let us find the corresponding Lagrangian L . Define the angular velocities from the special 
unitary matrices S: 

V a = iTk{Sd s &\ a }. (80) 

Then: 

L = ^V a 2 + lT,VdHl (81) 

Z a, z d 
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This Lagrangian gives the Hamiltonian in terms of the canonical momenta J a = J^-. There 
is a subtlety: upon quantization we have the generators of left and right SU(iV) rotations, 
related by the adjoint representation. The question is then to which the J's do correspond. 
This question is only relevant for the linear term,since the quadratic term is invariant under 
the adjoint. A short calculation shows that the left generators correspond to the J's. 

Let us now prove the following relation between the characters of the group and the inte- 
grated propagator of the SU(iV) symmetric top: 



J dVL (nu f \ exp {i(s f - Si )H} \toUi) 



= ^Trp^C/t)} exp[i(s f - Si )E (R,l)]. (82) 

R,l 

On the r.h.s. we integrate over SU(iV) matrices Q, with the measure normalized to 1. 
To prove this, insert the set of intermediate states \R; l,r > from Eq. [TTJinto the l.h.s. : 



dtt (ttU f \ exp {i{s f - Sl )H} \nUi) 

dtt Yl d R (toUf\R; I, r) exp [i(s f - Si)E°(R, I)] (R; I, r\QUi). (83) 

R,l,r 

The right index r can be summed over, and using Eq. [J_7J one finds: 

Y,{nU f \R-l,r){R-l,r\mJi) = D R (pU f (pUi)%. (84) 

r 

Integration over Q gives, for fixed I, using Eq. [7J3 

/ dQ D R {nU f {nU i )\ l = -^Tr {D R {U f U})} (85) 

J CLr 

Plug this result into the r.h.s. of Eq. IHS1 and we have the main result Eq. IK21 

Now we want to project out from the main result the representation R Q . This is done by 

letting / — > in the energy exponent and determining the minimum of E°(R,m). One such 

minimum is realized by = md and C 2 (R) = C 2 (Ro). 

So the sum over all irreducible representations R in Eq. IH21 reduces in this limit to the 

representation Rq singled out by the its corresponding highest weight H°. 

Strictly speaking there is only one minimum when H is the weight of the fundamental 

representation, or of any of the representations corresponding to a fully antisymmetric Young 

tableau. The reason is that for fixed A/"-ality C^-R) takes its minimum value for R being fully 

antisymmetric. 

For higher representations, like e.g. in SU(2) with weight jo( the number of boxes in the 
Young tableau), there is the representation R = jo — 1, with m = jo — 1, that minimizes E° 
as well. This is why in the original work [IH] the asymmetric top was taken, to provide an 
independent inertia in front of the Zeeman splitting in Eq.[7ji]and taking them independently to 
zero. So from now on the irrep Rq stands for one of the totally antisymmetric representations. 
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Of course the result contains the rapidly varying phase factor 



exp 



l{Sf - Si 



E (Ro, H 



21 



cxp 



i(sf- Si ) C2OR0) 



\^{ H l} ) I- 1 



= F(R 



According to appendix A the value in this minimum is (^(-Ro) — |Tr {Hq} = I^T 1 wi(N — 1)IN 
with W[ the weights defined by the Young diagram of R . This same factor F(R ) appears also 
in front of the path integral transcription of the matrix element: 



{flU f \ exp{i(s f - Si )H} \QUi) = F(R Q ) / DS(s) exp -i / dsL 



(86) 



with L as in Eq. |^ This formula, together with the fact that the left generators in the matrix 
element are correctly represented by the path integral can be proved straightforwardly. We 
introduce the periodic fluctuation variable Q(s) with Q(sij) = Q and transform the fluctuation 
variable S in Eq. EE! 

S{s) = tt(s)U(s) with DVt = DS. (87) 
The path integral now becomes: 

— i I dsL 



Dtt(s) exp -% J dsL (88) 

with L as in Eq. |HU but with the substitution Eq. |H7| in the angular velocities V a in (|5U|) : 

V a = Tr {Q(s)U{s)d s (Q(s)U{s)y X a } 

= Tr {(Q{s){U{s)d s U{ S y)Q{sy + Q{s)d s Q{ S y) X a } . (89) 
Finally use Eq. EHlto write the angular velocity as the gauge transformed potential A s : 

V a = Tr {{Q(s)(~tgA s )Q( S y + Q(s)d s Q( S y) X a } = Tr {(fi( S )V s fi(s) t ) A a } . (90) 

Remember V s fi(s) f = g .0 - igA)n(s) j . 

The final form of the Wilson loop in the irrep Rq and with highest weight Hq follows from 
Eq.|HUIB21and| 



w, 



lim 

l-*o 



DQ exp 



i 1 - j ds Vl - jds Tr Ih (VV ■ ^fi f ) 



(91) 



The line integral on the r.h.s. can be easily transformed into a surface integral, because it 
now involves the gauge transform average of potential projected on the highest weight, i.e. an 
Abelian potential. The surface S bounded by the loop L is covered by a set of nested loops 
L(r) with L(r = 0) = L and L(r = 1) shrunk to a point. Then, in an obvious notation: 



- j ds Ti{H (ttV s &)} 
J dr d r j> ds Tr {H nV s ^} 

J dr j> ds Tr {tf (e rs V r fiV s n f ) + e rs fiV r .V s fi f ) } 



(92) 



which gives then Eq. using the commutator [V r ,V s ] = —igF r s . The gauge transform in 
the last expression is now extended to all of the surface S and beyond. This surface is of course 
arbitrary, apart from its boundary. 
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Appendix C 

We give the explicit form of our Polyakov loops with jV-ality k. Let P(x, y, t) = ]ln=i U x {x + 
na,y,t), where U x (x) stands either for the original link variable U x (x) or a fuzzy |2SJ version 
of it with the same gauge transformation properties. Then our operators are 



(*)( f ) = a J2 (fc) (x, y + ma, t) (93) 



where the operator is one of 



k = 2: ^Tr{P 2 } ^(Tr{P}) 2 (94) 

k = 3: ^Tr{P 3 } ^(Tr{P}) 3 ^Tr{P 2 }Tr{P} (95) 

^ = 4: ^Tr{P 4 } ^(Tr{P}) 4 ^(Tr{P 2 }) 2 (96) 
^Tr {P 3 }Tr {P} ^Tr {P 2 }(Tr {P}) 2 , 

where the argument of P is the same as that of in Eq. EB1 The correlation functions we 
measure are (O^(0)(O^(t))*). 
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l:(3 = 115.0 
V = 24 x 28 x 24 


k = 1 


k = 2 


fc = 3 


k = 4 


L = 24 


1.536(19) 


2.67(13) 


/ 


/ 



= 115.0 
V = 16 x 20 x 24 


fc = 1 


fc = 2 


fc = 3 


fc = 4 


L = 16 
L = 20 


1.0075(31)* 
1.2787(33)* 


1.729(30) 
2.174(59) 


2.216(66) 
2.55(23) 


2.32(12) 
3.04(33) 




IIL/3 = 115.0 
y = 12 x 16 x 24 


fc = 1 


k = 2 


fc = 3 


fc = 4 


L = 12 
L = 16 


0.7338(38)* 
0.9980(45)* 


1.275(24) 
1.702(28) 


1.680(27) 
2.157(52) 


1.789(36) 
2.423(93) 



IV: (3 = 138.0 


fc = 1 


fc = 2 


fc = 3 


fc = 4 


V = 20 x 24 x 24 










L = 20 


0.808(10) 


1.436(15) 


1.826(20) 


1.985(29) 


L = 24 


0.991(13) 


1.745(20) 


2.214(42) 


2.28(11) 



V:0 = 138.0 


fc= 1 


fc = 2 


fc = 3 


fc = 4 


V = 16 x 20 x 24 










L = 16 


0.6348(58) 


1.1580(96) 


1.464(12) 


1.540(15) 


L = 20 


0.8253(75) 


1.462(15) 


1.845(19) 


1.937(29) 



Vl:p = 172.5 
V = 24 x 28 x 36 


fc= 1 


fc = 2 


fc = 3 


fc = 4 


L = 24 
L = 28 


0.5869(32) 
0.6961(58) 


1.0260(82) 
1.242(12) 


1.303(18) 
1.564(34) 


1.401(19) 
1.576(57) 




VII:/? = 172.5 
V = 20 x 24 x 36 


fc= 1 


fc = 2 


fc = 3 


fc = 4 


L = 20 
L = 24 


0.4883(31) 
0.5878(43) 


0.8939(65) 
1.042(16) 


1.109(12) 
1.308(23) 


1.213(18) 
1.381(22) 



Table 1: The masses of flux-tubes of different jV-alities. Values followed by a * were extracted 
from simulations employing a one-level algorithm. In addition, a /3 — 138 run on a 32 x 32 x 36 
lattice was done with am k=1 = 1.3347(49) and am k=2 = 2.326(49). 
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(3 = 115.0 


(3 = 138.0 


P = 172.5 




0.2558(6) 


0.2059(6) 


0.1582(12) 


7i 


0.624(78) 


0.67(12) 


0.27(19) 


02/01 L=00 


1.690(52) 


1.709(37) 


1.695(45) 


0-3/0-1 \ L=OQ 


2.01(12) 


2.162(72) 


2.175(88) 


04/01 L =00 


2.28(20) 


2.31(11) 


2.08(12) 



Table 2: Top: the fundamental string tension and the string correction coefficient (to be 
compared with the bosonic string value 75 = 71/6 ~ 0.5236. Bottom: the ratios of fc-torelon 
masses, extrapolated to L = 00 assuming 1/L 2 corrections at three different lattice spacings. 
The continuum limit a — > of these ratios are given in section 16.2.31 



(3 = 115.0 


k = 2 


k = 3 


A; = 4 


L= 12 
£=16 
L = 20 
L = 24 


1.694(33) 
1.688(28) 
1.686(46) 
1.728(84) 


2.213(36) 
2.141(51) 
1.98(18) 

/ 


2.352(49) 
2.319(91) 
2.35(26) 

/ 


Mean 


1.692(28) 


2.184(36) 


2.349(49) 



(3 = 138.0 


fc = 2 


jfe = 3 


k = 4 


L = 16 
L = 20 
L = 24 
L = 32 


1.782(22) 
1.749(24) 
1.744(30) 
1.734(37) 


2.239(28) 
2.208(31) 
2.207(51) 


2.352(32) 
2.356(42) 
2.27(11) 


Mean 


1.744(24)* 


2.222(28) 


2.349(32) 




/3 = 172.5 


k = 2 


k = 3 


jfe = 4 


L = 20 
L = 24 
L = 28 


1.786(18) 
1.732(17) 
1.763(23) 


2.202(28) 
2.177(33) 
2.195(28) 


2.404(44) 
2.317(35) 
2.230(84) 


Mean 


1.743(17)* 


2.195(28) 


2.304(35)* 



Table 3: The effective ratios of fc-string tensions, corrected for finite- length effects assum- 
ing; Eq. (JEU with bosonic string coefficient; at three different lattice spacings. The values at 
different L are fitted by a constant to give the 'mean' value. The mean values were obtained 
by averaging values at all L, except for those appearing with an asterisk *, where the shortest 
torelon was dropped. The continuum limit a — > of these ratios are given in section 16.2.31 
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m fc=2 


iK=2 + m* k=2 ) 


2m k=1 


II: L — 16 

L = 20 


2.122(60) 
2.49(25) 


1.932(33) 
2.33(13) 


2.015(6) 
2.557(6) 


III: L = 12 

L = 16 


1.677(37) 
2.244(75) 


1.476(23) 
1.973(42) 


1.468(4) 
1.996(9) 


IV: L = 20 
L = 24 


1.886(24) 
2.245(59) 


1.661(17) 
1.995(31) 


1.62(2) 
1.981(3) 


V: L = 16 
L = 20 


1.380(39) 
1.72(10) 


1.272(23) 
1.598(54) 


1.27(1) 
1.65(2) 


VI: L = 24 
L = 28 


1.328(15) 
1.672(31) 


1.181(9) 
1.457(16) 


1.174(6) 
1.39(1) 


VII: L = 20 

L = 24 


1.109(28) 
1.296(22) 


1.002(14) 
1.174(10) 


0.977(6) 
1.176(8) 



Table 4: The first-excited k = 2 torelon mass. The roman numbers refer to the different runs 
whose parameters are given in Tab. ^ 
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60 



60 



5.55 



5.25 







3D SU(8) spatial torelon lineic mass (k=l) 




0.05 0.1 0.15 0.2 0.25 0.3 

l/g 4 L 2 

3D SU(8) spatial torelon lineic mass (k=2) 




0.05 



0.1 



0.15 

/g 4 



0.2 



1/R 4 L 2 



Figure 4: The lineic mass of spatial torelons in the k — 1 and k 
intercept on the vertical axis yields the string tension. 
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2 sectors at (3 = 138. The 
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3D SU(8) k=2 to k=l string tension ratio 




2.6 
2.5 
2.4 
2.3 
2.2 
2.1 
2 



0.02 0.04 0.06 0.08 0.1 0.12 
l/(°iL 2 ) 



3D SU(8) k=3 to k=l string tension ratio 




0.02 0.04 0.06 0.08 
l/^jL 2 ) 

3D SU(8) k=4 to k=l string tension ratio 




3=115.0 
3=138.0 
3=172.5 



0.12 



0.02 0.04 0.06 0.08 0.1 0.12 



l/(°iL 2 ) 

Figure 5: Final analysis: the extrapolation of the string tension ratios to infinite string length 
(only the filled points are included in the extrapolation). The result of the extrapolation is 
shown beyond the vertical axis. 
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Figure 7: Top: the binding energy of /c-strings per unit length, in units of <j\ and rescaled by 
a factor N, as a function of N. Bottom: the string tension ratio <7fc ^ iV/2 , rescaled by a factor 
jr. The lattice data for SU(4) and SU(6) is taken from|2o]. 
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